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Abstract 

Asymptotics of solutions of a perfect fluid when coupled with a cosmological 
constant in four-dimensional spacetime with toroidal symmetry are studied. 
In particular, it is found that the problem of self-similar solutions of the first 
kind for a fluid with the equation of state, p = kp, can be reduced to solving 
a master equation of the form, 

2F(q,k)^-G(q,k) q '(0 = t 

For k = and k = —1/3 the general solutions are obtained and their main 
local and global properties are studied in detail. 



PACS numbers: 04.20.Dw, 04.20.Jb, 97.60.Lf 
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I. INTRODUCTION 



Recently, self-similar solutions of the Einstein field equations have attracted much atten- 
tion, not only because they can be studied analytically through simplification of the problem, 
but also because of their relevance in astrophysics [1] and critical phenomena in gravitational 
collapse [2-4]. In the latter case, as in statistical mechanics [5,6], the self-similar solutions 
represent the asymptotics of critical collapse in the intermediate regions. 

Recently, self-similar solutions of a scalar field were studied and all such solutions were 
found analytically [7]. It was shown that some of the solutions can be interpreted as repre- 
senting gravitational collapse of the scalar field. During the collapse, trapped surfaces never 
develop, and as a result, no black hole is formed. This is consistent with the general theorem 
presented in [8]. Although the collapse always ends in spacetime singularities, it was found 
that these singularities are spacelike and not naked [9]. 

In this paper, we shall study self-similar solutions of a perfect fluid with toroidal sym- 
metry and the equation of state, p = kp, where A; is a constant, and p and p denote, 
respectively, the pressure and energy density of the fluid. In particular, in Sec. II we write 
the corresponding Einstein field equations and define apparent horizons and black holes in 
such spacetimes, while in Sec. Ill we present the dimensional analysis of the problem. In 
Sec. IV, following Barenblatt [5], we study the self-similar solutions of the first kind, and 
find general solutions for k = and k = —1/3. For other cases the problem reduces to 
solving a master equation of the form, 

2F(9,A04-G(?,W = 7- (i- 1 ) 

In Sec. V, we study the local and global properties of the solutions found in Sec. IV, while 
in Sec. VI, we present our main conclusions. 
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II. THE EINSTEIN-FLUID EQUATIONS 



The general metric for a four-dimensional spacetime with toroidal symmetry can be cast 
in the form, 

ds 2 = g ab (x°, x 1 ) dx a dx b - L V^ '* 1 ) (dO 2 + dip 2 ) , (a, b = 0, 1), (2.1) 

where x a G (— oo, oo) and all x a have the dimension of length, L. The coordinates x A = 
{9, ip}, (A = 2, 3) are dimensionless with the surfaces x A = and x A = 2n being identified. 
The constant L has the dimension of length. The spacetimes usually have three Killing 
vectors, given, respectively, by 

Clearly, the metric is invariant under the coordinate transformation, 

x° = x° (V°, x' 1 ) , x 1 = x 1 (V°, x' 1 ) . (2.3) 

The energy-momentum tensor (EMT) for a perfect fluid, , takes the form 

Tya, = (p + p) u^Uy - pg^ u , (2.4) 

where w M , p and p denote, respectively, the four-velocity, energy density and pressure of the 
fluid. Using gauge freedom (2.3) we choose the coordinates such that 

g 01 (x°,x 1 )=0, ^ = ((?oo) 1/2 <?. (2.5) 

Denoting x° and x 1 by t and z, respectively, we find that the metric can be written as 

ds 2 = e x ^dt 2 - e^dz 2 - Ve""**'*) (d6 2 + dp 2 ) . (2.6) 

In this note we consider asymptotics of solutions of a perfect fluid. Using the remaining 
(trivial) coordinate transformation t' — t + t Q , where to is a constant, we assume that a 
spacetime singularity always starts to form at t — (if there is any). Then, it can be shown 
that the Einstein-fluid equations, 



R»v - ^9»uR + = (p + p) u^u v - pg^ u , (2.7) 
can be cast in the form, 

2fJ>,tz ~ V>,t {H,z + \z) - V, z v,t = 0, (2.8) 
e x ' u {2 (X :ZZ - n^ zz ) + /i iZ (ii :Z + \ z + u jZ ) + \ z (A i2 - vj)} 

-2 {y tU + n tU ) - fi tt (2/i t + u t + A, t ) + u t (u t - \ t ) = 0, (2.9) 

P = \{e' v [^, zz - fi, x (3fi, z + 2u z )] + e" V,t (^t - 2* t ) } + A, (2.10) 

V = X - {e' x [An iU - n, t (3^, t + 2A, t )] + e~ v ^ (fi, z - 2\ z )} - A, (2.11) 

where A denotes the cosmological constant, which can be parameterized as 

A = ±i (2.12) 

To study the formation of apparent horizons and black holes in the spacetimes described 
by the metric (2.6), we first give the definition of such notions. Following [8,10] let us first 
introduce two null coordinates u and v via the relations 

du = J(t, z) (e x/2 dt - e u/2 dz) , dv = K(t, z) (e x/2 dt + e u,2 dz) , (2.13) 

where J(t, z) and K(t, z) satisfy the integrability conditions, 

d 2 u d 2 u d 2 v d 2 v 
dtdz dzdV dtdz dzdt 

Without loss of generality, we assume that they are all strictly positive, 



(2.14) 



J(t,z)>0, K(t,z)>0. (2.15) 
Then, it is easy to show that, in terms of u and v, the metric (2.6) takes the form 

ds 2 = 2e 2a(u > v) dudv - TZ 2 (u, v) (d9 2 + dip 2 ) , (2.16) 

where 

a(u,v) = --\n(2JK), K(u, v) = L e~ M/2 . (2.17) 



It should be noted that the metric (2.16) is invariant under the transformations 

u = u{u), v = v(v). (2.18) 

Using this gauge freedom, we assume that the metric has no coordinate singularities in the 
coordinates, u, v and x A . This, in particular, implies that a is finite except at some points 
or on some surfaces where the spacetime is singular. 
Introducing two null vectors, / M and n^, by, 

which are orthogonal to the two-surfaces S = {x x : t,z — Constant}, we find that the 
expansions of the null rays u = Constant and v = Constant are given, respectively, by 

e n = n^sT = 2e" 2CT ^. (2.20) 

Definitions. A two-surface, S, of constant u and v (or constant t and z) is trapped, 
marginally trapped, or untrapped, according to whether Q x Q n > 0, Q$ n = 0, or 9[9 n < 0. 
Assuming that on the marginally trapped surfaces we have 6i\ s = 0, an apparent horizon 
is defined as a two-surface H foliated by marginally trapped surfaces, on which 9 n \ H ^ 0. 
The apparent horizon is outer, degenerate, or inner, according to whether C u 9i\ H < 0, 
£u@i\h = 0' or £uQi\h > 0) future if 9 n \ H < 0, and past if 9 n \ H > 0, where C v (C u ) denotes 
the Lie derivative along the vector / M (n M ). We define a black hole by the existence of a 
future outer or degenerate apparent horizon [8,11]. 

It is interesting to note that the above definitions can also be given in terms of the vector 
K x , 

K X K' X = 2e- 2a K v K u = iftV<m. (2.21) 

Thus, a two-surface, S, of constant t and z is trapped, marginally trapped, or untrapped, 
according to whether lZ t \ is timelike, null, or spacelike. 



From Eq.(2.13) we find that 



= 1 -a/2 dt_ = l_ v/2 
du 2J ' dv 2K 

^ = --e-»/\ % = —e^l\ (2.22) 
du 2J ' dv 2K K J 



from which we have 



e n = 2e- 2 ^ = ^ (e'^n t - e -"/X) . (2.23) 



III. DIMENSIONAL ANALYSIS 

Since t and z all have dimensions of length, and A, v and u are dimensionless, we have 

[t] = [z] = L, [A] = [i/] = M = 1. (3.1) 

On the other hand, because the Ricci tensor has the dimension of L~ 2 , then from Eq.(2.7) 
we find that 

[p] = [p] = L-\ [l] = L. (3.2) 
Thus, from t, z and I we construct two dimensionless quantities, 

f—f (3.3) 

Since the metric coefficients are functions of t, z and Z, and A, v and are dimensionless, 
we must have 

A (t, 2, Z) = A (£, 77) , 1/ (t, 2, /) = 1/ (£, 77) , n(t,z,l) = fi(£,ri) , 

p(t,z,l) = ±*(£,Ti), p(t,z,l) = -^n(t,r)). (3.4) 

For any given function f(£,r)), we have 
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f,t = ~ (t,f,£-vf,v) > 

f,z = — 

/,« = ^/,«- ( 3 - 5 ) 

Inserting these expressions into Eqs.(2.8)-(2.11), we find that 

(S^ - V"*) = 0, (3.6) 
2 [£ 2 + - 2r lt (/i,r, 5 + ^) + V 2 (», m + Km) + 2 £ + v a)_ 
-2 - [f (/^ + \s) - ^ (/^ + \^)] 

+ [f - \a) - v (Kv ~ \r,)\ [f iy* - - ^ - /^)] 

+e A -" [2 - \#) - p 4 (A* + ^) - A, € (A* - ^)] = 0, (3.7) 
$ = e^ [4/i^-/i, e (3/i, f + 2^)] 

+e- A - [£ - 2^) - 77 (/^ - 2^)] ± V, (3.8) 
n = e~ x {4 (e 2 /i,« " 277^,^ + V 2 Hm + 

- ~ W,v) [f (3/x, ? + 2A iS ) - r/ (3/i >r , + 2\ <ri )]} 

+e-"lM 4 - 2A, € ) ± 4^, (3.9) 

where the + (-) sign corresponds to A > (A < 0). 

Note that in the present case we have four equations, Eqs.(3.6) - (3.9), and five unknowns: 
A, v, (/,, $ and II. Thus, to determine these five functions uniquely, we need to have one 
more equation, which is usually provided by the equation of state of the perfect fluid. In 
this paper, we consider the case 

p = kp, (3.10) 

where A; is a constant. 



IV. SELF-SIMILAR SOLUTIONS OF THE FIRST KIND 



In this section, we consider the asymptotic behavior of the previous equations as rj — > 0. 
According to Barenblatt [5] (See also [6]), self-similar solutions of the first kind are defined 
as the existence of the limits, 

VO = !im A(^, rj), i/„(f ) = lim u(£, rj), 

rj— >0 ?7 — >0 

/i (0 = lim ^,77), $ (0 = lim$(e,^,p), 
n (0 = lining 77). (4.1) 
Substituting these expressions into Eqs.(3.6)-(3.9) and Eq.(3.10), we find that 

2 W - W (A' + 1/ + //) = 0, (4.2) 
2£ 2 (/ + O + $ (// + 1/) K " 2//' - A') + 4] 

+e w [2 (//' - A") - !/ (// - A') - A' (// + A')] = 0, (4.3) 
4 (CV)' - [(3 + *y + 2 (A' - fci/)] 

-e A ^ {4V - /i' [(1 + 3fc)^' - 2 (A' - fci/)]} = 0, (4.4) 

$ = e - " [4/i" - // (3/i' + 2i/)] + e" A f V (// - 2i/) , (4.5) 

IT = e" A [4 (£V)' - eV (3/i' + 2 A')] + e"V - 2 A') , (4.6) 

where a prime denotes ordinary differentiation with respect to £. In writing Eqs.(4.2)- 
(4.6), we have omitted all zero subscripts for simplicity, however, the functions //, v and A 
appearing in Eqs.(4.2)-(4.6) are functions of £ only. 
Integrating Eq.(4.2), we obtain 

// = 1^+^)/^ (4.7) 

where Co is a non-zero and otherwise arbitrary constant. Inserting it into Eq.(4.4), we find 
that 

2e ->+A)/2 

co(l + fc) (£ 2 e--e A ) 

x {fV [2 + (1 + fc)^'] + e A [2* - (1 + fc)£A']} , (4.8) 
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for k ^ -1. When k = -1, substituting Eq.(4.7) into Eq.(4.4) and considering Eq.(4.7) 
yields 

where /i is an arbitrary function of £ only. It can be shown that for such solutions Eq.(4.3) 
is satisfied automatically, while Eqs.(4.5) and (4.6) yield 

$ = = n, (fc = -l), (4.10) 

that is, the corresponding spacetime is vacuum. In fact, it is not only vacuum but also flat, 
as the corresponding Riemann tensor vanishes identically. Thus, no self-similar toroidal 
solutions of the first kind exist with the equation of state p = —p. In the following we 
consider only the case where k ^ —1. 

The combination of Eqs.(4.2), (4.3), (4.7) and (4.8) yields, 

£ [e (2e A - £ V) e " - e 2X ] A' 

-A£ 4 e 2 " [3(l + k)£i/ + Z\' + 4] 

+ke 2X [(1 + k)iv' - (5 + 4A;)£A' + A(2k + 1)] 

+2k£ 2 e u+x [(1 + k)£v' + (2k + 3)f A' - 4ife] = 0, (4.11) 
£ 4 e 2iy {(1 + kfi 2 (2u" - 3u' 2 - u'\') + (1 + k)i [(5k - 8)u' - A'] 
+4(A;-2)} 

-?e v+x {(1 + A;) 2 £ 2 (2i/' + 2A" + v' 2 + A' 2 - I0i/X') 
+2(1 + fc)£ [(Ilk + 4)i/ - 7A'] + 8(4A; + 1)} 

+e 2A {(1 + k) 2 i 2 (2A" - 3A /2 - v'X') + (1 + fc)£ [ku 1 + (16fc + 3) A'] 

-4k(4k + 1)} = 0, (M-l). (4.12) 

From the above, we see that the problem of solving the Einstein field equations now reduces 
to solving Eqs.(4.11) and (4.12) for the functions of v and A. Once they are found, Eqs.(4.7), 



(4.5) and (4.6) will give the functions //, $ and II. Eqs.(4.11) and (4.12) are non-linear, and 
finding the solutions of such equations is, in general, very complicated. To start with, let us 
consider the case k — 0. 



A. k = 



When k — 0, Eq.(4.11) reduces to 



A' (£V - e A ) = 0, 



(4.13) 



which has two solutions, 



(i)A = A , (ii) A = z/ + 21n(0, 



(4.14) 



where A is a constant. However, it can be shown that in the latter case the corresponding 
spacetime is vacuum, <E> = = II. Thus, in the following, we consider only the first case, for 
which we can set A = without loss of generality. Then, Eq.(4.12) reduces to 



2£ 2 (£V - l) v" - e (See" + l) v' 2 - 8 (£V + l) + 1) = 0. 



(4.15) 



To solve this equation, we introduce the function y(£) via the relation 



y = £V-i, 



(4.16) 



for which Eq.(4.15) takes the form, 



y 



5y + 4 / = _2 
y> 2y(l + y) y ' f 



It can be shown that the general solution is 



2Y 



Y 2 - 1 



+ In 



Y -I 



Y + l 



+ b, 



where a and b are two integration constants, and 



(4.17) 



(4.18) 



Y = ie v ' 2 . 
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(4.19) 



In terms of Y, Eqs.(4.5) and (4.8) yield 

(F 2 - if 



$ = 2a 



3-=(^) 2 (F 2 -l) 2 , (* = 0). (4.20) 



B. k^O, -1 

When A = v + 21n(£), spacetime is vacuum for any value of k. Then, for k ^ 0, —1 we 
consider the following ansatz, 

A(0=K0 + 21n(0 + 9(0, ( 4 -21) 
where g(£) is an arbitrary function of £. Substituting into Eq.(4.11) we obtain 

v+w+k-^y+frm =a (422) 

To solve this equation, let us first consider the case fc = —1/3. 

1. k = -1/3 
In this case Eq.(4.22) has the solution, 

q(0 = 90, (4-23) 
where go is a constant. Inserting Eqs.(4.21) and (4.23) into Eq.(4.12), we find 

i 2 v" - 2£V 2 - %u' - 12 = 0. (4.24) 
There are two particular solutions, 

(1) i/W = -21n(0, (2) i/f = -31n(0. (4.25) 

To find the general solution of Eq.(4.24) we set v = A(£) + and then Eq.(4.24) reduces 
to 
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A" - 2A' 2 - jA! = 0, 



(4.26) 



which has the general solution, 



A(0 = -\\n\A, - A o e 



(4.27) 



where A Q and A 1 are two integration constants. Thus, for k = —1/3 we have the following 
general solutions, 

\^) = qo -Un\A 1 -A o e , 

v(Z) = -\ln\? fa - M 2 ) , 

2A 1 e- q °l 2 i 

coiAi-A*?) 1 ' 2 ' 

12A 1 (1 - e qo ) 



e^ 2 = 



$ = -3n 



e^(A 1 -A e) 1/2 ' 



(4.28) 



2. -1/3 
In this case, Eq.(4.22) has the solution, 

"(0 = ~Y^k ln(0 " ihk [q + Ak ln (1 " eQ)] + u °> (429) 

where vq is a constant. Inserting the above expression into Eq.(4.12), we find 

V(q)q"-9(q)q' 2 -^h(q)q' = 0, (4.30) 

where 

/(g) = (3k 4 + 10fc 3 + 12k 2 + 6k + l) (l - e 3q ) 

+ (6A; 5 + 23A; 4 + 34A; 3 + 24A; 2 + 8A; + l) e 2q 

- (3k 5 + 16A; 4 + 32A; 3 + 30A; 2 + 13k + 2) e q , 
g(q) = k (7k 4 + 22k 3 + 24k 2 + 10k + l) e 3q 

- (lOk 5 + 37k 4 + 50k 3 + 28k 2 + Ak - l) e 2q 
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+ (3k 5 + 12k 4 + 20k 3 + 18k 2 + 9k + 2) e q 
-3(1 + A:) 2 (l - A: 2 ) , 
%) = k (3k 3 + 7k 2 + 5k + l) e 3q 

- {6k 4 + 17k 3 + 17A: 2 + 7k + l) e 2q 
+ (3A; 4 + 13A; 3 + 19A; 2 + Ilk + 2) e q 

- (3k 3 + 7k 2 + 5k + l) . (4.31) 

It can be shown that when q'(£) = the corresponding solutions represent a vacuum space. 
When q' ^ 0, Eq.(4.30) can be written as 

2F(q)[\n(q')]'-G(q)q' = ^ (4.32) 

where 

Solving Eq.(4.32) is not a trivial exercise. One may first try to find a particular solution of 
it, say, q s (0- Once such a solution is known, setting q = q s (£) + <?o(0> we find that Eq.(4.32) 
reduces to the following form for go(0' 

[\n(q o ')}'-H(q o ,q s )q o ' = 0. (4.34) 

Unfortunately, we have not yet been able to find such a particular solution. 

V. PHYSICAL AND GEOMETRICAL INTERPRETATIONS OF THE 

SELF-SIMILAR SOLUTIONS 

In this section, we study the local, as well as the global, properties of the self-similar 
solutions for k = and —1/3, obtained in the last section. Note that, although these 
solutions were found by taking the limit i] — > 0, in this section we extend them to any 
t G (-oo,0). 
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A. Self-Similar Solutions With k = 



These are the solutions given by Eqs.(4.14)-(4.20). Rescaling the coordinates t and z 
and using the conformal transformation g' = B 2 g^ u , where B is a constant, without loss 
of generality, we set 

c = ^a, Ay = 1, (5.1) 

for which the metric reads 

ds 2 = dt 2 - ^dz 2 - (d9 2 + V) • (5.2) 
Then, the corresponding energy density is given by 

'=5** = w I(n (5 ' 3) 



where 

Y - 1 



/(F) = -b(Y 2 -l) + 2Y+ (Y - 1) In 



F+l 



(5.4) 



a(Y 2 - 1) 

From the above expressions one can see that the spacetime is singular at t = and £ = —a/6, 
where 

£(Y) = ~, when Y = 0, ±oo. (5.6) 

The singularity at t = is always spacelike, while the nature of the singularity at £ = —a/6 
depends on the values of Y. The normal vector to the surface £ = —a/6 is given by 

d { a \ a 



Then, 



n 2 t2 f a 2 > 0, F = ±oo, 



-oo < 0, F = 0. 
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Therefore, when Y = the corresponding spacetime singularity at x — —a/b is timelike, 
and when Y = ±00, it is spacelike. From Eqs.(2.17) and (5.2) we obtain 



L 



Thus, we have 



Y 2 - 1' 



From Eq.(5.5) we also find that 

d£(Y) 4a 



(5.11) 



dY P{YY 

which shows that £(Y) is a monotonically increasing (a > 0) or decreasing (a < 0) function 
of Y, depending on the sign of the constant a. To have the energy density p non-negative, 
the solutions are restricted to the following regions for different a, 

(>0, a>0, 

Y£ = (5.12) 
[<0, a<0, 

as one can see clearly from Eq.(4.20). To further study the properties of the solutions, 
it is convenient to consider the following four cases separately: (a) a > 0, b > 0; (b) 
a > 0, b < 0; (c) a < 0, b > 0; and (d) a < 0, b < 0. 
Case (a) a > 0, b > 0: In this case we find that 

' —00, Y = ±00, 
J(y) = j&, F = 0, (5.13) 
,±2, Y = ±l, 
and 7(F) = has two real roots, Y±, with the properties 

Y + > 1, -1< y_ < 0, (5.14) 

as can be seen from Fig. I. 
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FIGURES 
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FIG. 1. The function I(Y) versus Y for a > and 6 > 0. 



From Eqs.(5.4) and (5.5) we find that 

' -a/b, Y = 0, ±oo, 

£00 = jo, r = ±i, (5.15) 

k oo, Y = Y±. 

Then, one can see that £00 mus t behave as that given by Fig. 2, from which, together with 
Eq.(5.12), we can see that the energy density is non- negative only in the following three 
regions, 

(i) Y G [-oo, -1], or £ G [-a/b, 0], 

(ii) Y G (y_, 0], or £ G ( oo, 0], 

(Hi) Y G [1, F+), or £ G [0, oo). (5.16) 







i 5 




-1/ 
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i/ 


Y + Y 






o / 








-a/b 
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FIG. 2. The function £(Y) versus y for a > and b > 0. 



In the region Y G [— oo, —1] or £ G [—a/6, 0], the spacetime is singular at £ = —a/b or 
F = — oo. From Eq.(5.8) we can see that this singularity is spacelike. In the (t, z)-plane, 
this region is between the two lines £ = —a/b and £ = 0(z = 0), as shown in Fig. 3, in 
which it is referred to as Region /. The metric is singular at £ = or F = — 1. However, 
this singularity is a coordinate one, as one can see from the expression of p, which is finite 
there. The nature of this surface is null. In fact, introducing the normal vector N\ to this 
surface by 

dz 



dx x 



— °A> 



we find that 



N X N X 



Y 2 



o, 



(5.17) 



(5.18) 



as £ — > and Y — > —1. On the other hand, from Eq.(5.10) we find that in this region 1Z : \ 
is always timelike, that is, the whole region is trapped. 

Region //, where Y G (Y_, 0] or £ G (— oo, 0], is the region in between the two spacetime 
singularities t — and £ = —a/b. In this region, TZ t \ is always spacelike, that is, this region 
is untrapped. However, because of the two spacetime singularities, it is difficult to interpret 
this region physically. 

In Region III, where Y G [1, Y+), or £ G [0, oo), the vector 7£ A is always timelike, that 
is, all of this region is trapped. The spacetime is singular at t — 0. 





, t 


Singularities 




/ : 












i 


in 


(Q ! e n >o) 


(e,e n >o) 
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FIG. 3. The (t, z)-plane for a > and b > 0. The spacetime is singular on the lines £ = —a/6 
and £ = 0. 



Case (b) a > 0, 6 < 0: In this case we find 

' oo, Y — ±oo, 

/(F) = -|6|, y = o, 

,±2, Y = ±l, 
for which /(V) = has also two real roots, Y±, but now with 



o < y + < i, r_ < -l, 



as shown in Fig. 4. 





ioo i 


I / 






2 








-1 


/ 1 


Y 
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FIG. 4. The function /(y) versus y for a > and b < 0. 
Then, it can be shown that 



' a /\b\, Y = 0, ±00, 
0, F = ±l, 



00, Y = Y±, 



(5.19) 



(5.20) 



(5.21) 



and the curve of £(Y) versus Y is given by Fig. 5, from which we can see that the energy 
density is non-negative only in the regions, 



(1) Y E (Y_, -1], or £ G (-00, 0], 
(ii) Y G [1, 00], or £G [0, a/\b\}. 



(5.22) 
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FIG. 5. The function £(Y) versus Y for a > and 6 < 0. 



In the region Y E (Y_, —1] or £ G (— oo, 0], the spacetime is singular at t = or Y = Y_. 
In the (t, z)-plane, this region is the one where t < and z < 0, marked as Region I in Fig. 
6. The metric is singular at £ = or z = 0. As in the previous case, this singularity is a 
coordinate one, and the surface is null, as shown by Eq.(5.17). Moreover, in this region lZ t \ 
is always timelike, except for the line t — or Y — YL. Thus, in the present case the whole 
region is trapped. 

In Region II, where Y e [1, oo] or £ e [0, a/|6|], the spacetime is singular on the line 
£ = a/|&| or z = — at/\b\, and the nature of the singularity is spacelike, as can be seen from 
Eq.(5.8), considering the fact that now £ = a/\b\ corresponds to Y = oo. In addition, 1Z t x is 
also timelike, that is, Region II in the present case is also trapped. 



FIG. 6. The (t, z)-plane for a > and b < 0. The spacetime is singular on the lines £ = a/\b\ 
and t = 0. 




Singularities 



II 



(e,e n >o) 



(e,e n >o) 
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Case (c) a < 0, b > 0: In this case 7(F) is given by Fig. 1, and £(Y) behaves as that 
given by Fig. 7, from which we can see that the energy density is non-negative only in the 
regions, 

(i) Y e [-oo, -1], or £ E [\a\/b, 0], 
(it) Y E (Y_, 0], or £e [|a|/&, oo), 

(m)Ye[l, y+), or fe(-oo, 0]. (5.23) 













V 


lal/b 






Y_ 




Y + Y 


-1 ^ 











FIG. 7. The function £(Y) versus Y for a < and 6 > 0. 

The corresponding three regions of Eq.(5.23) in the (t, z)-plane are shown in Fig. 8. In 
particular, the spacetime is singular on the lines t — and £ = \a\/b, and Regions 7 and 777 
are trapped, while Region 77 is not. The metric is singular on the line z = or F = ±1, 
which is a null surface. 









Singularities 















III 






i 


(e, e n >o> 




0, 


e„>0) 



FIG. 8. The (t, z)-plane for a < and b > 0. The spacetime is singular on the lines £ = |a|/6 
and i = 0. 
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Case (d) a < 0, b < 0: In this case I(Y) is given by Fig. 4, and £(Y) behaves as that 
given in Fig. 9, from which we can see that the energy density is non-negative only in the 
regions, 

(i) Y G (Y_, -1], or f G [0, oo), 

(ii) y G [0, y+), or £g (-oo, -|a/6|], 

(m) y G [1, oo], or [-|a/6|, 0]. (5.24) 





j 








\-l 


Y+ 


\l Y 
V 


Y_ 


\ 








-la/bl 







FIG. 9. The function £(Y) versus Y for a < and b < 0. 

The corresponding three regions of Eq.(5.24) in the (t, z)-plane are shown in Fig. 10, 
where the spacetime is also singular on the lines t — and £ = —\a/b\, and Regions I and 
III are trapped, while Region // is not. The metric is singular on the line z = 0ory = ±l, 
and this line is null. 





t 


Singularities 




/ 










in 


i 


(e, e n >o> 


(6,e n >0) 



FIG. 10. The (t, z)-plane for a < and b < 0. The spacetime is singular on the lines £ = \a/b\ 
and t = 0. 
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B. Self-Similar Solutions With k = 



-1/3 



These are the solutions given by Eq.(4.28). Without loss of generality, we set 

A 1 = 1, c = 2e~ qo/2 . (5.25) 
Then, the solutions are given by 

A = g - r In (l - a£ 2 ) , 



2 

v — — In 
2 

" = ln (T^i)- < 526) 



«" (i - , 



for which we have 



'=-*=?(r^P' (5 - 27) 

where p = (1 — e 90 )/(4e 90 ) and a = A . Thus, when a > 0, the solutions are valid only in 
the region where |£| < £o as shown in Fig. 11, and the spacetime is singular at t — and 
at £ = ±£ = ia^ 1 / 2 . These singularities are all spacelike. This can be seen easily when 
t — 0. For the cases £ = ±£ , we first introduce the normal vectors by 

Then, we find that 



n±n±g^ = e - (u+qo) - e qo j ~ e-^KoJ+w) (i _ e *>) > o. (5.29) 

Thus, the singularities on £ = ±£ are indeed spacelike. On the other hand, it can be shown 
that now we have 



n = Lo- 



rn 



U a K' a = ^- -r, (1 - e qo ) > 0. (5.30) 

t 2 e(a/e-i) 1/2 

Thus, the whole region — x < £ < x is trapped [cf. Fig. 11]. It should be noted that the 
metric is singular on z = or £ = 0, but it can be shown that it is a coordinate one. 
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Singularities 




FIG. 11. The spacetime in the (t, z)-plane for a > 0, for which the solutions are valid only in 
the region |£| < £o- The spacetime is singular on the lines £ = ±£o, which are always spacelike. 
For a < 0, the solutions are valid in the whole half plane t < 0, and the spacetime is singular only 
on the line t = 0. 



When a < 0, the solutions are valid in the whole half plane t < 0, and the spacetime is 
singular only on t — 0. It can also be shown that the whole region t < is now trapped, as 
one can show that 6i9 n > for all time t < 0. 



VI. CONCLUSIONS 

In this paper, we have studied the asymptotics of solutions of a perfect fluid when coupled 
with a cosmological constant in four- dimensional spacetimes with toroidal symmetry. We 
found that the problem for self-similar solutions of the first kind with the equation of state, 
p = kp, can be reduced to solving a master equation of the form, 

2F( q ,k)^-G(q,k)M) = j. 

Although we were not able to solve this equation for all k, we did obtain the general solutions 
for k = and k = —1/3. The local and global properties of these solutions were studied in 
detail, and it was found that no apparent horizons develop during the evolution of the fluid, 
although trapped regions indeed exist. This is consistent with the general theorem obtained 
previously [8]. 
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Finally we note that spacetimes with toroidal symmetry are locally indistinguishable 
from those with plane symmetry. In fact, by first unwrapping the angular coordinates 9 
and ip and then extending the ranges to x A e (—00, 00), we obtain spacetimes with plane 
symmetry. With a different identification of coordinates, one may obtain topologies other 
than toroidal and plane geometries. 
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